Atomic Squeezing under Collective Emission 

V.I. Yukalovi'2'4 and E.P. Yukalova^'^ 



^Bogoluhov Laboratory of Theoretical Physics, 
Joint Institute for Nuclear Research, Dubna 141980, Russia 

^ Nonlinear Optics Division, Institute of Physics, 
Adam Mickiewicz University of Poznan, Poznan 6I6I4, Poland 

^Department of Computational Physics, Laboratory of Information Technologies, 
Joint Institute for Nuclear Research, Dubna 141980, Russia 

^Institut fiir Theoretische Physik, 
Freie Universitdt Berlin, Arnimallee I4, D-14195 Berlin, Germany 



Abstract 

Atomic squeezing is studied for the case of large systems of radiating atoms, when 
collective effects are well developed. All temporal stages are analyzed, starting with the 
quantum stage of spontaneous emission, passing through the coherent stage of super- 
radiant emission, and going to the relaxation stage ending with stationary solutions. A 
method of governing the temporal behaviour of the squeezing factor is suggested. The 
influence of a squeezed effective vacuum on the characteristics of collective emission is 
also investigated. 



PACS numbers: 42.50.Fx, 42.50.Dv 
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1 Introduction 



Quantum engineering of squeezed states has been a subject of great interest during the past 
several years. The squeezed states of hght [1-3] and atomic squeezing [4-6] have attracted 
the major attention, though the effect of squeezing can be defined for any two operators 
[7]. Atomic squeezing is directly related to the radiation field squeezing [8] and, vice versa, 
squeezed atomic states can be created by irradiating atoms with squeezed light [9-12]. Res- 
onant finite-level atoms are usually represented by spin operators, because of which atomic 
squeezing is often termed spin squeezing. Since the spin operators, employed for describing 
finite-level atoms, do not represent real spins, but rather are convenient mathematical tools, 
one also uses for atomic squeezing the names of dipole squeezing [1,4] or pseudospin squeezing 
[5]. Recently this type of pseudospin squeezing has been studied for Bose-condensed atoms 
with two internal states [13], for two-component mixtures [14], and for multimode Bose- 
Einstein condensates [15-17]. There is a variety of possible practical applications of atomic 
squeezing, e.g., for atomic spectroscopy and atomic clocks [6], for atom interferometers [18], 
and for quantum information processing and quantum computation [19]. 

Atomic squeezing for radiating atoms is usually studied by considering a single atom 
placed in the so-called squeezed vacuum, that is, subject to the action of squeezed light. 
Then one considers the influence of squeezed vacuum on the radiation process of the atom, 
whose relaxation characteristics can be noticeably modified by the presence of the vacuum 
squeezing [20,21]. An effective squeezed vacuum can be generated, e.g., by parametric am- 
plifiers [20,21]. Squeezing of light can be achieved under its interaction with nonlinear media 
[22] , such as Kerr medium [23] . When light interacts with condensed matter, there can also 
arise squeezed exciton-polariton and phonon-polariton states [24-26] and squeezed magnon 
states [27]. Squeezed states of transverse acoustic phonons interacting with hght were ob- 
served in experiment [28]. From the theoretical point of view, squeezed polariton states 
appear under the canonical Bogolubov transformation diagonalizing the matter-radiation 
Hamiltonian [24,29]. 

If not a single, but several atoms are considered, then one has to take into account their 
effective interaction owing to the photon exchange through the common radiation field. This 
interaction may also infiuence atomic squeezing. The problem of resonance fluorescence from 
two atoms driven by a squeezed vacuum field was analysed in Rcfs. [30-34]. 

The aim of the present paper is to consider atomic squeezing for an opposite case, when 
the number of radiating atoms is large. Then, contrary to the case of one or two atoms, 
collective effects come into play. Owing to photon exchange, strong interatomic correlations 
may arise. We aim at studying the intrinsic atomic squeezing in a multiatomic radiating 
system. It is especially interesting to analyze the temporal behaviour of squeezing, starting 
from an initially inverted nonequilibrium system, which begins with spontaneous radiation, 
passing through the incoherent quantum stage to the stage of well developed coherence, and 
then going to the relaxation stage. The evolution of atomic squeezing through all these 
stages is the primary concern of the present paper. Understanding the temporal behaviour 
of atomic squeezing, we may find a way for its regulation. One such a way for governing the 
evolution of squeezing is suggested. 

Another problem, we consider here, is the influence of squeezed vacuum on the radiation 
characteristics of a multiatomic system. Such a setup, as is mentioned above, can be prepared 
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by incorporating the radiating atoms into a medium with squeezed polariton states. We show 
that squeezed vacuum can substantially change the system attenuation. 

Thus, we study both sides of the story, the temporal evolution of atomic squeezing in the 
presence of collective correlations and the influence on this evolution of squeezed vacuum. 



2 Atomic squeezing 

In this brief section, we give the main definitions, introducing those quantities that will be 
investigated in what follows. 

The notion of squeezing can be defined for any two operators, say A and B, which are 
not necessarily self-conjugate. The definition stems from the Heisenberg uncertainty relation 

A^A) A^B) >^\< [A,B] >f , 

in which the dispersion of an operator is given by 

A\A) = < A+A > -\<A>W 

with A"*" being the Hermitian conjugate of A. One says that A is squeezed with respect to 
B if 

A\A)<l\<[AB]>\. 
This naturally suggests to introduce the squeezing factor 

2A^(A) 

quantifying how A is squeezed with respect to B. If so, then Q{A, B) < 1. In terms of factor 
(1), the Heisenberg relation reads as 

Q{A,B) Q{B,A) > 1 . 

Resonant atoms are characterized by means of spin operators. For an i-th atom, these are 
the ladder operators Sf", representing atomic transitions, and the operator S'f , representing 
the population difference. Here and in what follows, the operators S'j^ = S'"(rj,t) are treated 
as the Heisenberg operators associated with spin 1/2, which corresponds to two-level atoms. 
The Heisenberg equations of motion for these operators translate into the evolution equations 
for the related statistical averages. We define the transition functions 

Ui = 2<Sr>, u; = 2<S+>, (2) 

local coherence intensity 

Wi = A\<Sf>W (3) 

and the population difference 

Si = 2<St> . (4) 
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The introduced quantities are the functions of space and time, Ui — u{ri,t), wi — w{ri,t), 
and Si = s(rj, t). 

The operators Sf and arc the most convenient for characterizing atomic squeezing. 
Deahng separately with the operators S^, S^, and ^'f , one would meet the so-called trivial 
squeezing caused by rotation [35]. For the considered operators, one has 



AW) = ^(1±«.)- \ ^\St)^ \ (l-.^) , \<St>\ = \ 
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The squeezing factor, not involving trivial squeezing [35], and introduced according to defi- 
nition (1), writes as 

Q{St,Sf)^^-^. (5) 



One says that the atomic state is squeezed if Q < 1. Prom the physical point of view, this 
means that, under Q < 1, the population difference can be measured with a higher accuracy 
than the transition characteristics, such as the coherence intensity (3). To investigate the 
behaviour of the squeezing factor (5), one has to know functions (3) and (4). 

3 Collective emission 

Consider a system of N resonant atoms inside matter. The appropriate Hamiltonian reads 
as 

H^H, + Hf + Haf + H^f . (6) 
The Hamiltonian of resonant atoms, with a transition frequency cuq, is 

N .1 



i=l 



Ha^T.^o[^ + S!) . (7) 



Here and in what follows, we employ the Gaussian system of units and set the Planck 
constant h = 1. In the field Hamiltonian 

E is electric field and H = V x A is magnetic field. Por the vector potential A, the Coulomb 
gauge 

V ■ A = (9) 
will be used. The atom-field interaction is represented by the dipole Hamiltonian 

^ /I \ 
^a/ = -E -J^-A^ + Pi-E° , (10) 

i=l ^ 

in which the transition current is 

Ji^iu;o{dSt-d*S-) (11) 
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and the transition dipole operator is 



Pi = dSt + d*S- , (12) 

with d being a transition dipole. Here Aj = A(rj,t) is the vector-potential operator and 
= Eo(rj, t) is an external electric field. The last term in Eq. (6) represents the matter-field 
interaction described by the Hamiltonian 

Hmf \ j Jmat(r, t) ■ A(r, t) dr , (13) 

where jmat(r,t) is a local density of current in matter. 

Writing down the Heisenberg equations of motion, we shall use the standard commutation 
relations 

■£;"(r,f), Af'{r\t)]^Amc5af3{r-r') , 
in which c is light velocity and 

2 1 
Sa(3{r) = g Saf3 S{r) - — D^p{r) 

is the so-called transverse delta- function [22], where 

= ;3 ' 

with = r"/r and r = |r|, is the dipolar tensor. 

The common way of treating the Heisenberg equations of motion is to take their statistical 
averages, with the semiclassical approximation, which results in the system of the Maxwell- 
Bloch equations [36]. However, for the present consideration, it is more useful to follow 
another approach, by eliminating the field variables and obtaining the equations containing 
only the spin variables [37,38]. This method is, first, more convenient, since for calculating 
the squeezing factor (5), we need exactly the averages of the spin variables. And, second, this 
approach allows for the consideration of the quantum stage of collective atomic emission, 
while the Maxwell-Bloch equations are classical, thus being unable to treat the quantum 
stage of evolution. 

From the Heisenberg equations of motion for the field variables, using the Coulomb gauge 
(9), we have 



where the density of current is 



r(r,t) = E |E'^a/3(r - r.)jf (t) + / 5«Mr - r')jL,(r',t) dr' 

U=i J 



(15) 



The general solution to Eq. (14) writes as 



I r — r' I \ dr' 



(16) 
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where A^^c is the solution to the related uniform equation, which represents vacuum fluctu- 
ations. 

Counting time from the initial moment t = 0, we keep in mind that all processes are 
defined for t > 0, which can be properly taken into consideration by the causal condition 

snt)^o {t<o). 

Also, one should not forget that the very notion of resonant atoms, having a well-defined 
transition frequency, has sense only when the interaction of radiation with an atom is sub- 
stantially weaker than intra-atomic interactions. This makes well grounded the Born ap- 
proximation 

ST (^t- ^) = Sr{t) Q{ct - r) 6""°' , (17) 

important for treating the retardation effects in Eq. (16). Here ko = cuq/c and ©(•) is a 
unit-step function. 

Taking into account the atomic self-action yields, as is known [39], to the appearance of 
the natural width 

70 = ^ \d\'kl . (18) 

This is usually generalized by including in the equations of motion the longitudinal, 71 = 
1/Ti, and transverse, 72 = I/T2, relaxation parameters. 
The vector potential (16) can be written as a sum 

A = A„(jc -|- A.^ad + ^dip + ^mat ■ (19) 

Here A^^c is the vacuum vector potential. The long-range part of the vector potential 
generated by radiating atoms is 

A.^(M)=i:3^j,(*- ^) ■ (^o) 

The short-range part is due to the dipolar term resulting in 

AyM)^-i:i:/^^^/(-^)-'. (21) 

And the last term in Eq. (19) describes the local vector potential produced by matter, 

A^„,(r, ^) = E / ^^^TT^ fmat (r", t - ^) dv' dv" . (22) 

All short-range parts can be combined into the variable 

{(r, t) =2kod- {A^ac + Adip + A^at) ■ (23) 

There are in the system two different spatial scales. One is the radiation wavelength 
A = 2ttc/ujo and another, the mean interatomic distance a. At optical frequencies, one 
usually has a <C A. The long-range potential (20) acts at the distance of A, while the typical 
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distance of action for the short-range variable (23) is a. Because of the existence of two 
different scales, we may employ the scale separation approach [38,40-42]. Then we separate 
the operator variables in two parts. One part is responsible for long-range interactions related 
to the long-range potential (20). Since the latter is expressed through the spin operators, 
we may call the set 

{S} = {S^,S-,St\i^l,2,...,N} 
the family of long-range variables. And the set of short-range terms 

{0^{e(r,t)|reV}, 

where V is the system volume, can be referred to as the family of short-range variables, or 
short-range fluctuating fields. Treating {S} and {^} as two different operator families, we 
may introduce two types of averaging for any function F = F{{S}, {^}) of these variables. 
To this end, we define the statistical averaging over the spin degrees of freedom as 

< F > = TV{5} P F{{S}, {C}) , (24) 

with p being a statistical operator. While the averaging over the fluctuating flelds is denoted 
as 

« F » = Tr{^y p F{{S}, {C}) . (25) 

Then we write the Heisenberg equations of motion for the spin operators S~ — S~{ri, t) 

and Sf = S^{ri,t) and average these equations over the spin degrees of freedom as dcflncd 
in Eq. (24), which results in the equations for functions (2), (3), and (4). It is convenient 
to pass from the summation to integration by means of the replacement 

where p is the atomic density. In this way, the equations for functions (2), (3), and (4) are 
reduced to the equations for the functions u — u{r, t), w — w{r, t), and s — s(r, t). To write 
down these equations in a compact form, we deflne the effective force 

/(r, t) = /o(r, t) + fradir, t) + ^(r, t) . (26) 

The flrst term here is due to the external held, 

/o(r,i) = -2id-Eo(r,i). (27) 

The second term is caused by the radiation potential (20), 

fradir, t) = 2ko < d • Aradir, t) > . (28) 

And the last term is the fluctuating fleld (23). The radiation force (28), taking into consid- 
eration Eq. (20), can be presented in the form 

frad(r,t) = -i^oP I icir - r',t) u{r',t) - ^ G*(r - r',t) u*(r',t)\ dv' , (29) 
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with the transfer function 

exp(zfcor) 
kor 

As a result, we obtain the equations 

— = -{iuQ + 72) + /s , — = -272W + (li*/ + s , 
(9s 1 

^ = -2 («7 + Ai)-7iG^-C), (30) 

where is a stationary population difference for a single atom. When there is no external 
nonresonant pumping, then ^ = — 1. Generally, allowing for pumping, one has ( in the 
interval —1<C<1. Equations (30) are stochastic differential equations because of the 
fluctuating field ^ entering the effective force (26). 

The evolution equations (30) can be further simplified in the frame of the scale separation 
approach [38,40-42]. For this purpose, we notice that there are different temporal scales, 
which follows from the existence of the small parameters 

(31) 

UJq UJq COq 

with the assumption that the atom-field interactions are much weaker than Uq. Then Eqs. 
(30) tell us that the function u varies in time much faster than w and s. Hence, the slow func- 
tions w and s are temporal quasi-invariants with respect to the fast function u. Therefore, 
we can invoke the averaging technique [43] generahzed to the case of stochastic differential 
equations [38,41,42]. 

For concreteness, let us consider a sample having typical of lasers cylindrical shape. The 
axis of the cylinder is along the 2;-axis, which is selected by the propagating seed field 

Eo(r, t) = - El e*('=^-'^*) + - E* e'^^*^^-'"*^ , (32) 

in which k = cu/c, and which is in resonance with the atomic transition frequency, 

lAI 



< 1 , (A = u; - a;o) ■ (33) 



In a realistic physical situation, the wavelength A = 2'nc/uj is much smaller than the radius, 
and length, L, of the cylindrical cavity, 

^ « 1 , ^ « 1 ■ (34) 

When considering a small-aperture cylindric sample, whose Fresnel number F = /XL 
is less than one, then there exists the sole propagating mode and the single-mode description 
becomes valid. If the Fresnel number is larger than one, then there arise several transverse 
modes, represented by photon filaments [44,45]. The effective Fresnel number for each fila- 
ment is less than one. Therefore the problem again can be reduced to the consideration of 
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separate filaments, for each of which the single-mode description is applicable [46-49]. In 
the single-mode picture, we have 



u{r, t) = u(t) e'"^ , w(r, t) = w{t) , s(r, t) = s{t) . (35) 

We substitute Eqs. (35) into Eq. (30), multiply the latter by e~*'^^, and average over 
space. We also take into account that the interaction time Tint = a/c is very small, since for 
a ~ 10~^ cm, one has Tint ~ lO"-*^^ s. Therefore after the very short dynamic interaction stage 
< t < Tint long-range interactions are already established. Being so short, this dynamic 
stage can be neglected, and one may look for the solutions for t > Tint- 
In order to simplify the following expressions, let us introduce the effective coupling 
parameters 

7o f sin(/cor - kz) , 70 /" cos(/cor - kz) 

g = p— dr , g = p— dr , 36 

72 J kor 72 J kor 

where the integration is over the sample volume V — ttR^L. These parameters characterize 
the effective strength of interaction between the system of atoms and radiation field. We 
also introduce the collective width 



collective frequency 

and the dynamic detuning 



r = j,{l-gs), (37) 
fl = LLio + g'j2S , (38) 



S = u-Q = A- g'-f2S . (39) 
Using the above notation, we can represent the solution for the fast variable as 



u 



( - e-(^^+r)* + e-*-* + s f i{t') e-(^^+r)(*-0 dt' , (40) 

V 6-\-iVJ d-\-iV Jo 



which is obtained from the first of Eqs. (30), keeping fixed all temporal quasi-invariants, 
and where 

,,,d.E, (41) 

and the effective fluctuating field is 

J C{r,t) e-^'^ dv . (42) 

Then we substitute solution (40) into the second and third of Eqs. (30) and average them 
over time and over the fluctuating fields, keeping in mind that the latter are zero-centered, 
such that -C 3>= 0. We define the effective attenuation 

^ '''^''^ (l-e-^*)+r3, (43) 



52 + r2 
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in which the first term is caused by the seed field (32), while the last term is the quantum 
attenuation 

Ta = Re ^lim ^ dt < C{tW) > e-(^"+^)(*-*') dt' (44) 

due to the action of the fiuctuating field ^{t), playing the role of an effective vacuum. As a 
result, we obtain the evolution equations for the slow variables 

^ ^ -2^^(l - gs) w + 2f s\ ^ = -^72«;-fs-7i(s-C) . (45) 

These are the final equations defining the functions w and s that we need to know in order 
to calculate the temporal behaviour of the squeezing factor (5). 



4 Stages of evolution 

Equations (45) are rather general and allow one to consider various situations. Here we 
concentrate on the case that is, to our mind, the most interesting, when radiation coherence 
develops in the system as a self-organized process. This means that at the initial time no 
coherence is imposed upon the atomic system by external forces, but atoms are assumed to 
be inverted, which implies the initial conditions 

w(0) = , s(0) = 1 . (46) 

In order to study the intrinsic atomic squeezing, happening in the radiating system 
without a special preparation of a squeezed vacuum, we center, first, our investigation on 
the white-noise vacuum, when the fluctuating flelds are correlated as 

«r(i)e(0» = 2735(^-0- (47) 

With Eq. (47), we get from Eq. (44) the equality T3 = 73. 

Aiming at studying a self-organized process, we should not include external fields acting 
on atoms. Recall that the role of the seed field (32) has been to select a resonant atomic 
mode and to prescribe the cylindric symmetry of the system. After this has been done, we 
may set ui 0. Consequently, Eq. (43) yields F = Fa = 73. 

As can be checked by straightforward calculations [46-49], for A ^ a, one has g ^ 1, 
where the coupling g is defined in Eqs. (36). Keeping this in mind, we shall simplify the 
following formulas by assuming large g ^ 1. 

Consider now the relaxation process for the system of inverted atoms, when the initial 
conditions (46) are satisfied. Recall that, as is mentioned above, there always exist a very 
short interaction stage, when all atoms radiate independently, not yet feeling each other. 
This stage lasts in the time interval 

< t < Tint {interaction stage) , (48) 

where the interaction time Tint — o/c. Being so short, the interaction stage can be neglected. 
The following dynamics is described by Eqs. (45). 
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In the next stage, atoms begin feeling each other through the photon exchange of the 
common radiation field. However atomic phases are not yet correlated, because of which 
radiation is incoherent. This is the quantum stage governed by local quantum fluctuations. 
The stage lasts in the interval 



till the crossover time 



{quantum stage) 



To 



tc = — ■ 

2y 



At this stage, the solutions of Eqs. (45), with F = 73, are 

w ~ 2-fst , s ~ 1 - (71 + 73 - 7iC)t . 



(49) 
(50) 

(51) 



Coherence in the system develops in a gradual way, being negligible before the crossover 
time (50), but becoming important after this time. After the time tc, the coherent stage 
comes into play. Equations (45), in the interval tc <t <^Ti, yield 



w — sech^ 
where the delay time is 

and the pulse time is 



s — tanh I 

9 



Try 



tn = 



— 'P 



1 + ln 

= 1l 
— 

9 



Ag 



72 
73 



(52) 

(53) 
(54) 



The coherence intensity w{t) reaches its maximum at the delay time (53), when 

w{to) = 1 , s{to) = ^ . 

Coherent radiation happens in the time interval between tc and to + ^p, after which the 
coherence intensity w{t) exponentially decays. Hence the coherent stage is characterized by 
the interval 

tc < t < to + Tp {coherent stage) . (55) 
The decay of coherence after this stage is described by the solutions 



w ~ 4 exp 



+ 

9 



(56) 



where to t -C Ti. This signifies the beginning of the relaxation stage, which can be defined 
as occurring at times 

to + Tp < t < Ti {relaxation stage) . (57) 

When there is no external stationary pumping, so that ( = ~1, the relaxation stage is rather 
trivial. But if there is a strong pumping, such that C ~ 1) the relaxation stage consists of a 
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series of coherence bursts with gradually diminishing amplitudes. The analysis of Eqs. (45) 
shows that the coherent flashing is not periodic though it is possible to find the approximate 
time between neighbouring coherent bursts which is 



The number of bursts can be estimated as 

--^(^)''- 

Finally, for very large time, the process ends with the quasistationary stage, when 

Ti <^t {quasistationary stage) . (60) 

Defining the stationary solutions of Eqs. (45), we keep in mind that, by definition, < w < 1 
and — 1 < s < 1. Of possible fixed points, we select those that satisfy these restrictions. 
Generally, when there are several fixed points, it is necessary to accomplish the Lyapunov 
stabihty analysis and to choose the stable fixed point. Following this way, we find the stable 
stationary solutions to Eqs. (45) as 

w* = ^^^^ {iM - 1) + isib - 73) + 71 [^^ - 73(2 + ^C)]} , 

= ^^(^ + 90 + ls-b 
271^ 

in which 

62 = [7i(l + K)+73]'-47iV- 

These solutions essentially depend on the value of the pumping parameter C. There are three 
qualitatively different situations: no pumping or weak pumping (C f=i — 1), medium pumping 
(C ~ 0), and strong pumping {C ^ 1). 

For the case of the medium pumping (C ~ 0), from Eqs. (61), to third order in C, we get 

^* _ 7i73C^ ^ 7i73(7i - 73)gC^ 
(71 + 73)^72 (71 + 73)^72 

s* ~ 7iC _ 7i73^C^ _ 7?73(7i - 73)g^C^ /g2) 
~ 7i + 73 (71 + 73)^ (71 + 73)^ 

When ICI Ki 1, we may simplify Eqs. (61) taking into consideration that g ^ 1. Then to 
second order in 1/g, we find 

* _ (71 + 73)|C| + (71 -73)C _ / X (71 +73)|Cl + (71 -73)C 

272^ ^^'^^'^ 27i725^|C| 

. _ C- ICI , (7i + 73)|CI + (71 -73)C _ 73 ( . 

2 + 27i^|C| 7i^^lCI ■ ^ ^ 
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For weak pumping, when C < 0, this yields 



^.^73^ _ (71 + 73)73 ^ + (64) 



And for strong pumping, when ^ > 0, we have 

* 7iC 71 +73 * 1 73 

w — , s ~ — . (65) 

129 129^ 9 li9 Q 



These solutions allow us to calculate the behaviour of the squeezing factor at different stages 
of the radiation process. 



5 Dynamics of squeezing 

Atomic squeezing is characterized by the squeezing factor (5), which can be written as 

1 - 

Q = ■ (66) 



If the atomic system is prepared in the inverted state, with initial conditions (46), at the 
interaction stage (48), one has 

g ~ (i ^ 0) . (67) 

This means a very strong atomic squeezing. 

At the quantum stage (49), using solutions (51), we have 

g ^ ^i^i-^l±^ x/2t . (68) 



There is yet atomic squeezing, which diminishes with time, since Q grows with time. 

Coherent effects, arising at the coherent stage (55), destroy squeezing. Thus, t — Iq, 
the squeezing factor can be approximated as 

g ~ 1 - 4 = ^0) , (69) 

9 

from where g ~ 1, since g ^ 1. 

At the relaxation stage (57), the squeezing factor, first, grows as 

g ~ - e*/^^ , (70) 
9 

and then either monotonely tends to a stationary value, if there is no pumping, or experiences 
a series of oscillations, when the pumping is present. 

At the quasistationary stage (60), the squeezing factor approaches its stationary value 
depending on a particular physical situation. Thus, for the case of medium pumping (^ a; 0), 
employing solutions (62), we find 



g^ 



71 + 73 _ (71 - 1^)9 

^ 71 Id 2(71 + 73), 



(71) 
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Since \(\ — > 0, the squeezing factor is large, hence there is no squeezing. 

When there is no pumping, or it is very weak, such that C ~ — 1, then using Eqs. (64), 
we get 

While if the pumping is strong, that is ~ 1, then from Eqs. (65), we find 

7i + 73 + (71 - 73)2 



Q 




1 + 



2715 ^lig 



2^2 



(73) 



Strong pumping destroys squeezing. If the pumping is absent, so that — —1, one has 
atomic squeezing for sufficiently large coupling parameters, when 

Q<1 (^>^)- (74) 

The temporal behaviour of the squeezing factor (66) can be determined more accurately 
by solving numerically Eqs. (45). For this purpose, we consider these equations for F = 73. 
To simplify notation, wc measure time in units of 7^^, and the relaxation parameters 71 and 
73 in units of 72. Then Eqs. (45) reduce to the form 

dvj ds 

— = -2(1 -t/s)w + 2735^ , = -^«^-73S-7i(s-C) ■ (75) 

Solving Eqs. (75) numerically, we calculate the squeezing factor (66). The results are shown 
in the following figures, in all of which the initial conditions (46) are assumed. 

When the stationary external pumping is absent, so that C, — —1, the temporal behaviour 
of the coherence intensity w{t) and population difference s{t) is qualitatively the same for 
different parameters 71, 73, and (?, although varying these parameters, of course, modifies 
the quantitative values of w and s. The typical dependence of w{t) and s{t) on time is 
presented in Fig. 1. The related squeezing factor Q{t) is shown in Fig. 2. The coherence 
intensity, proportional to the coherent atomic radiation intensity, exhibits the characteristic 
superradiant bursts at the delay time to, when the population difference falls fast from s — 1 
to almost s = — 1. The squeezing factor is a nonmonotonic function of time, whose shape 
essentially depends on the system parameters. It may possess a maximum and a minimum, 
as in Fig. 2, just one maximum, two maxima and a minimum, as in Fig. 3, or several 
maxima and minima. 

In the presence of an external pumping, when C = 1) the regime of pulsing superradiance 
evolves. Then the functions w{t) and s{t) display oscillations, as is shown in Fig. 4, which 
results in the appearance of a series of oscillations of the squeezing factor in the related Fig. 
5. The number of these oscillations and their sharpness can be varied by changing the system 
parameters. For instance, the number of oscillations, according to Eq. (59), increases with 
increasing the coupling g. This is demonstrated in Figs. 6 and 7, which should be compared 
with Figs. 4 and 5. 

Even more possibilities in regulating the temporal behaviour of the squeezing factor arise 
in the regime of punctuated superradiance, when the self-organized evolution of the radiating 
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system is punctuated by short vr-pulses inverting atoms to their initial inverted state s = 1. 
The regime of punctuated spin superradiance was suggested in Ref. 50. The situation for 
atomic systems is analogous to that for spin assemblies. Realizing the regime of punctuated 
superradiance for the considered atomic system, we can radically alter the behaviour of the 
squeezing factor as is demonstrated in Figs. 8 and 9. 

In this way, atomic squeezing in the system of radiating atoms can be regulated by 
means of different techniques. One possibility is by preparing a system with the parameters 
assuring the desired behaviour of the squeezing factor. Another method is by organizing the 
regime of punctuated superradiance, with a chosen temporal location of inverting pulses. 
The combination of these methods is also admissible. 

6 Squeezed vacuum 

As is discussed in the Introduction, when resonant atoms are inserted into a medium with well 
developed polariton effects, then this medium can form an effective squeezed vacuum [22-28]. 
Gardiner [20] noticed that spontaneous emission of a single atom can be essentially distorted 
by the influence of squeezed vacuum. Here we study how squeezed vacuum influences the 
attenuation of an ensemble of many radiating atoms. 

An effective squeezed vacuum appears due to the interaction of electromagnetic field with 
matter. Local fluctuations of this effective squeezed vacuum are superimposed on the local 
fluctuations caused by dipolar atomic interactions. The latter are well represented by white 
noise. The effective squeezed vacuum is usually modelled by parametric oscillators [20-23] . 
Therefore the ffuctuating fleld (42) can be written as a sum of two terms: 

m - m + E 7K) (&.e— + 6;e-^*) . (76) 

Here the flrst term corresponds to white noise, for which 

« Ut) » = , « CoitMt') » = 2735(t - t') , (77) 

analogously to Eq. (47). The second term in Eq. (76) describes a system of oscillators 
modelling an effective squeezed vacuum. The creation, and annihilation, bg, operators 
satisfy the Bose commutation relations. The spectrum cuq — cu-q, in which q is momentum, is 
positive. The spectral function 7(0;) is assumed to be symmetric with respect to the central 
line cus, so that 

The averages of the oscillator operators give 

«:65» = 0, <&5&p> = V9P' <6q6j> = (l+n5)5qp, (78) 

where ng is a momentum distribution function. The averages ^ hqhp ^ and ^ fejfej^ ^ in 
the case of normal vacuum are zero, while for a squeezed vacuum such anomalous averages 
do not vanish, 

< h\hl > = mg/^{ujg + ujp- 2ujs) , (79) 
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where A{uj) is the discrete delta-function 

and rriq is prescribed by the properties of squeezed vacuum. The momentum distribution Ug 
is semipositive, but the anomalous factor niq is, in general, complex, 

Uq = n{ujq) > , mq = \m{ujq)\e''^'' , (80) 

with the phase (fg defined by the squeezing field proportional to cos{2u!st + (ps)- The white- 
noise and squeezed-vacuum fluctuations are not correlated with each other, implying that 

< Ut)bg > = . 

Characteristic expressions for Ug and rriq are given in Appendix. 
For the correlator of the fluctuating field (76), we find 

« emit') » = 2735(i - t') + Ef(^.) [n,e-''(*-*')+ 

+(1 + nje-^'^''^*-*') + ^^e^'^''*+^(^'"^ + ^^e-^'"^*-^^^'"'-'"'')*'] . (81) 
Substituting this in Eq. (44) , we obtain the quantum attenuation as a sum 

Ta = 73 + r^es + ^non + ■ (82) 

Here 73 is due the white noise caused by dipolar interactions. The resonant part of the 
quantum attenuation, due to the effective vacuum fluctuations, is 



n[uj) 



(1 - e-^^*) , (83) 



(52 + r2 

where ui is the frequency of the seed field, F and 5 are defined in Eqs. (37) and (39), respec- 
tively, and the detuning from resonance is assumed to be small, \S\ -Co;. The nonresonant 
attenuation writes as 

Ug l + ria 



(84) 



{n-LOgy + T^ (n + LUg)^ + T^ 

And the last term in Eq. (82) is the attenuation generated by squeezing, 

T. I ^ M 2/ s Fcos<^s + 2a;,sin(^s _rt /orx 
F, = -|m(u;)|7 (u;) ^ ^2 ^ ■ (^5) 

At the initial time t = 0, the resonant attenuation (83), as is seen, vanishes. Hence, 
the start of the radiation process is mainly due to the noise attenuation 73, nonresonant 
attenuation (84) and the squeezing attenuation (85). For large times, the resonant attenua- 
tion (83) becomes more important than nonresonant quantities (84) and (85). In this way, 
the presence of squeezed vacuum results in the change of the attenuation 73 to an effective 
attenuation (82). Changing the attenuation leads to the shift of the delay time (53), defining 
the maximum of the first coherent burst, and also to the variation of the stationary solutions 
(61) to (65). Consequently, there would arise some modifications in the temporal behaviour 
of atomic squeezing characterized by the squeezing factor (66). The overall evolution of the 
latter remains qualitatively the same. But the details of this evolution can be regulated by 
the properties of squeezed vacuum. 
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7 Conclusion 



We have considered the atomic squeezing for an ensemble of many resonant atoms. The used 
approach makes it possible to analyze all stages of emission, which are: interaction stage 
(48), quantum stage (49), coherent stage (55), relaxation stage (57), and quasistationary 
stage (60). The temporal behaviour of the squeezing factor is studied in the course of its 
evolution through all these stages. Both analytical estimates and numerical calculations are 
presented. We show how it is possible to govern the behaviour of the squeezing factor in the 
regime of punctuated superradiance. This opens the way of an effective regulation of atomic 
squeezing by suppressing the fluctuations in the population difference with respect to the 
fluctuations of the transition dipoles. 

In the media with a well developed polariton effect, an effective squeezed vacuum can 
arise. By incorporating resonant atoms into such media, one may change the characteristics 
of the atomic system. The influence of the squeezed vacuum on the collective radiation 
by atoms consists in modifying the temporal characteristics of emission and varying the 
stationary solutions. Hence, by preparing different kinds of squeezed vacuum, it is feasible 
to regulate collective atomic emission and the evolution of atomic squeezing. 

Acknowledgement 

We are very grateful to R. Tanas for many useful discussions and helpful advice. Financial 
support from the Bogolubov-Infeld Program and the German Research Foundation (DFG 
grant Be 142/72-1) is appreciated. 



17 



Appendix 



To present explicit examples of the momentum distribution n{w) and of the anomalous 
term \m{u) \ in the case of an effective squeezed vacuum, we can resort to the standard case of 
a parametric oscillator [20-22], with a squeezing field of intensity e (in dimensionless units). 
One employs the notation 

l^ = l-e , i/ = 7 + £ , 

with 7 being a cavity damping rate. Then, for a nondegenerate parametric oscillator, one 
has 



n{u}) 



+ 



1 1 

+ 



(A, + + 1/2 

1 



+ 



+ 



1 



where A^ = ui—uig and the parameter n characterizes a two- mode squeezed field, representing 
the displacement from the central frequency of squeezing, where the two-mode squeezed 
vacuum is maximally squeezed. If the parametric oscillator is weakly nondegenerate, such 
that \k\ <^ |As|, or degenerate, with k — 0, then 



n{u!) 



|m(a;) 



A2+//2 Aj + u\ 

( 1 1 



A2 + /i^ A2 



As is evident, generally, the normal, n(a;), and the anomalous, |m(a;)|, terms are of the same 
order of magnitude. 
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Figure Captions 



Fig. 1. Temporal behaviour of the coherence intensity w (sohd hne) and population 
difference s (dashed line) for 71 = 0.1, 73 = 1 (in units of 72), and g = 100, without pumping, 
so that C = —1. Both w and s are dimensionless. Time is measured in units of 7^^. 

Fig. 2. Evolution of the squeezing factor Q for the parameters of Fig. 1. Time is 
measured in units of 7^^. 

Fig. 3. Squeezing factor Q for the parameters 71 = 73 = 1, 51 = 100, and ^ — —1. Time 
is measured in units of 7^^. 

Fig. 4 Oscillations of the coherence intensity w and population difference s as functions 
of time (in units of 7^^) in the presence of pumping, when ^ = 1, for the parameters 
71 = 73 = 1, 9' = 10. 

Fig. 5 Evolution of the squeezing factor Q for the parameters of Fig. 4. Time is 
measured in units of 7^^. 

Fig. 6. Regulating the number of oscillations of w and s, under ^ = 1, by varying the 
system parameters, which are 71 = 1, 73 = 0.1, and g — 100. Time is measured in units of 
72-^ 

Fig. 7 Evolution of the squeezing factor Q for the parameters of Fig. 6. Time is 
measured in units of 7^^. 

Fig. 8 Governing the behaviour of the squeezing factor Q in the regime of punctuated 
atomic superradiance for the parameters 71 = 0.1, 73 = 1, g = 100, and ( = —1, with 
equidistant inverting pulses at the moments of time ti — 0.1, t2 = 0.2, = 0.3, ^4 = 0.4, 
and ^5 = 0.5. Time is measured in units of 7^^. 

Fig. 9 Governing the evolution of the squeezing factor Q in the regime of punctuated 
atomic superradiance for the same parameters as in Fig. 8, but with nonequi distant inverting 
pulses at the moments of time ti — 0.1, ^2 = 0.2, ts — 0.3, ^4 = 0.4, = 0.5, = 1, = 1.1, 
^8 = 1-2, ^9 = 2, tio — 2.4, and tu — 2.5. Time is measured in units of 7^^. 
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